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We consider a mesoscale model for nano-sized metaparticles (MPs) composed of a central sphere decorated 
by polymer chains with laterally attached spherocylinder. The latter mimics the mesogenic (e.g., cyanobiphenyl) 
group. Molecular dynamics simulations of 100 MPs reveal the existence of two novel morphologies: uCol^ 
(hexagonal columnar arrangement of MPs with strong uniaxial order of mesogens collinear to the columns 
axis) and wCol^ [the same arrangement of MPs but with weak or no liquid crystalline (LC) order]. Collinearity of 
the LC director and the columnar axis in uCol^ morphology indicates its potentially different opto-mechanical 
response to an external perturbation as compared to the columnar phase for the terminally attached mesogens. 
Preliminary analysis of the structures of both phases is performed by studying the order parameters and by 
visualisation of the MPs arrangements. Different mechanisms for the mesogens reorientation are pointed out 
for the cases of their terminal and lateral attachment. 
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1. Motivation 

Colloid particles, polymers and LC molecules represent main building blocks of soft matter physics 
m When two or more of such blocks are combined into a MP, the latter forms a meta-material that ex¬ 
hibits a range of new morphologies and new effects that are not observed for any of its pure constituents. 
Examples to mention are: LC elastomers, LC dendrimers, decorated nanoparticles and others EHZl, all 
of these having already found a number of applications in thermo- and photo-controlled elasticity, plas- 
monic resonance, photonics and medicine. Most applications rely on the symmetry of the equilibrium 
morphology, in which each constituent is “responsible” for particular property of a meta-material. For 
instance, its elasticity is usually controlled by a polymer subsystem, optical properties are governed by 
both the behaviour of mesogens and the arrangement of gold nanoparticles (if any). Therefore, the type 
of mutual arrangement of the constituent parts is crucial regarding the new potential applications of each 
particular meta-material. 

Let us concentrate on the MPs built out of a spherical core decorated by polymer chains (spacers) 
each ending by a mesogen. The core mimics either a solid nanoparticle (e.g., gold nanoparticle (6][7l) or 
averaged in time shape of a dendritic scaffold Ellsj. It was found experimentally that the most important 
aspects of decoration are as follows: the surface density of mesogens on the outer shell of MP, the length 
of a spacer, and the exact way mesogens are attached chemically to the spacer ItHTTI. In particular, 
both terminal and lateral attachment can be realised chemically (H |5l IMTTl [depicted schematically in 
figure(a) and (b)] and the difference between the morphologies observed in these two cases are in the 
focus of this study. 
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The case of a terminal attachment of mesogens [figure (a)] has been studied more widely out of 
the two (Hzl. In particular, at low grafting density of chains such MPs adopt a rod-like conformation 
and self-assembly into lamellar smec¬ 
tic A morphology [see, figure (a)]. 

With an increase of a grafting den¬ 
sity, the disc-like conformation for MPs 
is favoured over the rod-like one, and 
the MPs self-assemble into columns 
which arrange themselves hexagonally 
[see, figure |^(b)]. Both structures pos¬ 
sess uniaxial symmetry denoted via the 
symmetry axis d. For the morphology 
(a) it is defined as a normal vector to 
the layers (and it also coincides with 
the preferential direction for MPs rods), 
whereas for its (b) counterpart d is de¬ 
fined along the columns of MPs. Despite 
this common feature, the different type of MPs arrangement in these morphologies will affect their elastic 
properties and their behaviour in plasmonic resonance applications (if the cores represent gold nanopar¬ 
ticles). Their behaviour will differ even more in the applications related to the optical activity and orienta¬ 
tional order of mesogens. Indeed, morphology (a) is characterised by an uniaxial preferential orientation 
of mesogens defined via nematic director n, whereas for the case of morphology (b) the global direc¬ 
tor does not exist. The distribution of mesogens orientations in this case is fiat radial centered around 
the axes of each column [as indicated in figure [^(b)]. One can clearly see some analogy with the case 
of LC elastomers, where the role of a symmetry axis d can be attributed to the preferential orienta¬ 
tion of the polymer backbones |2l|3l. It is well known that their opto-mechanical applications depend 
heavily on the mutual arrangement of d and n axes, the main- and side-chain architectures being an 
example l2l[3l [T2WT^ . 

Therefore, a question arises whether 
some other morphologies can be ob¬ 
tained for the MPs with lateral attach¬ 
ment of mesogens [figure (b)] which 
are characterised by different mutual 
arrangement of the d and fi axes, lead¬ 
ing potentially to new optical and opto¬ 
mechanical behaviour. Based on this 
MP geometry, one would expect the for¬ 
mation of hexagonal columnar mor¬ 
phology depicted in figure (c). It 
should be quite similar to its counter¬ 
part (b) (shown in the same figure) 
in terms of its mechanical properties 
and the arrangement of cores, but, cru¬ 
cially, is characterised by uniaxial ori¬ 
entational order with the nematic director shown as fi therein. Moreover, are there any other morpholo¬ 
gies possible? We try to provide the answer in the following section by means of coarse-grained molecular 
dynamics simulations. 



Figure 2. (Color online) (a) and (b): schematic representation of 
lamellar smectic and hexagonal columnar morphology, respec¬ 
tively, observed for the MP decorated by polymer chains with 
terminally attached mesogens. (c): uniaxial hexagonal colum¬ 
nar morphology that is expected to exist for the case of laterally 
attached mesogens. 
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Figure 1. (Color online) Schematic representation of a MP made 
of a spherical core (shown in pink) decorated via spacers 
(shown in gray) each ending by a mesogen (shown in blue), 
(a) MP with terminal attachment of mesogens, (b) the same 
with lateral attachment, (c) and (d): reorientation options for 
attached mesogens for both models (see text for more details). 


2. Results and discussion 

The model for a MP with laterally attached mesogens is based on its counterpart with terminal at¬ 
tachment described in detail in references (T7l[T8l. The central sphere represents a core (be it a gold 
nanoparticle or a dendritic core) decorated by grafting Nch = 32 chains to it. We use the annealing-like 
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grafting, when the first bead of each chain slides freely on the surface of a central sphere. The effective 
diameters of soft beads are based on the coarse-graining of the atomistic model for the generation 3 LC 
dendrimer performed in reference (19). These are: 2.14 nm for a central sphere, 0.62 nm for the first bead 
and 0.46 nm for the remaining beads of the spacer which connects mesogens to the central sphere. 

The spheres interact via soft repulsive potentials of a quadratic form (T7]|l8). The mesogens are mod¬ 
elled as soft spherocylinders of the breadth D = 0.374 nm and of the length-to-breadth ratio L/D - 3. 
These dimensions are also used for the visualization purpose. The MP is kept together by means of har¬ 
monic spring potentials, with the same parameters used for the model with terminal attachment I17II18I . 
As far as each monomer is assumed to represent three hydrocarbons 23 and their dimension is less than 
the Kuhn segment length, we introduce the pseudo-valent angles which account for the effective spacer 
stiffness I17II18I . 

Orthogonal lateral attachment of the mesogen to the spacer [figure[^(b)] is maintained by employing 
two potentials. First one is the usual harmonic bond potential between the center or the last monomer of 
the spacer and the center of the spherocylinder, the bond length is 0.3 nm. The second potential is also of 
a harmonic type with respect to the angle between the latter bond and the long axis of a mesogen 1^ . 
The equlibrium value for this angle is set equal to Oq = nl2. The force constant for the latter potential 
is set equal to its counterpart for the pseudo-valent angle potential mentioned above (the absolute value 
can be found in references I17II18I ). 

In this study, we consider the MPs with A/ch = 32 grafted chains only. The initial configuration is build 
out of 100 such MPs arranged randomly in a simulation box. The initial orientation of grafted chains is 
radial, with respect to the central sphere. Typical box dimensions are of the order of 15 20 nm and the 
periodic boundary conditions are used in the simulations. First, a short run is performed of 10000 molecu¬ 
lar dynamics steps with a small time step of At = 2 fs to eliminate initial interparticle overlaps. Following 
the scheme employed in references (17] [HI, we performed field-assisted self-assembly simulations first. 
These are done at a fixed temperature T = 520 K by imposing an external orienting field. The latter is 
introduced via additional energy term ^ cos^df, which contains a contribution from each 

ith mesogen. Here, di is the angle between the field direction and the orientation of ith mesogen, F is the 
magnitude of the field. Throughout this study we assume that the field is always directed along the Z axis 
of the Cartesian frame. Its magnitude F sets a time-scale for the mesogens reorientation and we found 
empirically the values in the range of F = (2-^4) • 10~^® J to be optimal for our state points. The simulations 
are performed in anisotropic isobaric ensemble, where the temperature and each of the three diagonal 
components of the stress tensor are maintained constant by means of the thermostat and three separate 
barostats 11311161 . This makes possible the auto-adjustment of the simulation box shape which at the end 
becomes commensurate with the pitch of the emerging ordered phase. One should remark that the po¬ 
tential interactions used in this model are predominantly soft repulsive (save for the mesogen-mesogen 
potential introduced by Lintuvuori and Wilson (21)). Therefore, to maintain the density required for the 
ordered phases, the system needs to be stabilized by the pressures of P ^ 50 atm I17II18I . 

The field-assisted self-assembly simulations are performed for the range of pressures P = 50-^200 atm, 
the duration of each run is 10® molecular dynamics steps with a time step of Af = 20 fs. For all values of 
the pressure within this interval, the system indeed self-assembles into the uniaxial hexagonal columnar 
morphology depicted schematically in figure[^(c) and thereafter referred to in this study as uColh- Depen¬ 
dence of this morphology on the pressure is refiected only in the density of stacking of MPs into columns, 
as well as in the total density of the system. We will provide simulation snapshots for this morphology 
below. The next stage is to test the stability of this structure and the changes it undergoes in a wide range 
of temperatures when the aiding field is removed = 0). To provide a quantitative analysis of these 
changes, let us introduce a set of the order parameters. 

The level of orientational order of mesogens is described by the nematic order parameter defined as: 

5 n = {P2icos , ( 1 ) 

where 6i is the angle between the orientation of ith mesogen and nematic director axis h, and P 2 ix) = 
I (3x^ -1) is the second Legendre polynomial. The averaging is performed over all mesogens in a system 
and over the time trajectory denoted as (.. The initial director was set along the Z axis in the course 
of a field-assisted run. And, similarly to the case of MPs with terminally attached mesogens (17]|18), we 
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found a negligible drift of the director away from this axis for all temperatures where the orientation 
order exists whatsoever. Therefore, in this study we assume 6i to be the angle between the orientation of 
ith mesogen and the Z axis. 

The amount of the “hexagonality” 
in the arrangement of MPs in the 
XY plane can be characterized by the 
hexagonal order parameter Su, intro¬ 
duced as follows. Let us consider ith 
MP and identify its first coordination 
circle in the XY plane. To this end, we 
evaluate the bond vectors Ifjc = = 

,/f,} between its center to the 
centers of each of its neighbours in¬ 
dexed via k. The bond lengths between ith and kth MPs projected onto the XY plane are: = 

neighbour k is assumed to belong to the first coordination circle of ith MP if the 
following conditions fulfill: Rmin < < ^max and |/^^| < Zmax- For our model, we choose the following 

set of parameters: Rmin = 4.7 nm, Rmax = 5.3 nm and Zmax = 0.2 nm but the choice, obviously, depends on 
the geometry of the model MP. After the first coordination circle that contains MPs are identified, the 
arbitrary axis in the XY plane is chosen (e.g., the X axis) and the angles are evaluated for each kth 
MPs belonging to this circle [see, figure(a)]. Then, both the local hexagonal order for ith MP and 
its global counterpart Sr can be defined via equation ^ as: 



Figure 3. (Color online) (a) Flat hexagonal cluster of MPs and 
the definition of bond angles within it. (b) A column of MPs 
stacked along Z axis. 
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Here, j = \/^. The same equation also contains a definition for the other order parameter, Sq, which 
accounts for the “columnarity” of the uColh phase, i.e., the amount of stacking of the MPs discs along 
the Z axis. It can be related to the average columns height which is identified by considering ith MP 
and drawing through its center an imaginary cylinder of radius Rq that extends along the Z axis [see, 
figure [^(b)]. Then, one takes into account the number of MPs so that their centers are found inside 
this cylinder (periodic boundary conditions are not taken into account). Normalisation factor A^nax is 
introduced for convenience, to bring the range of Sq values to the same interval as the other two order 
parameters, 5n and Sr. We used the following values: i?c = 1 nm and Nmax - 15. These notations explain 
the expressions for the columnar order parameter Sq in equation 
Let us concentrate on the behaviour 
of the order parameters S^, Sr and Sq 
when the uniaxial hexagonal columnar 
phase uColh, obtained by the field-aided 
self-assembly, is equilibrated at various 
temperatures. The plots demonstrating 
the changes are provided in figure]^ At 
relatively low temperatures, T < 350 K, 
the values of all three order parame¬ 
ters are essentially non-zero and one 
obtains the uniaxial hexagonal colum¬ 
nar phase uColh first depicted schemat¬ 
ically in figure [^(c). It is characterized 
by both uniaxial nematic order along 
the Z axis and regular columns of disc¬ 
shaped MPs aligned along the same axis 
and arranged hexagonally in the XY 
plane (see snapshots in the left-hand 
column of figurej^. With an increase of 
the temperature, all order parameters 



300 350 400 450 500 550 T, K 

Figure 4. Changes undergone by the order parameters Sr, Sr 
and Sc when uColh morphology, obtained by the field-aided 
self-assembly, is equilibrated at various temperatures from 
300 K to 560 K. The lines connect data points and serve as 
guides only. 
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gradually decay to their minimum values. First of all, let us note that 5 n decays very fast (practically 
linearly), in contrast to the case of terminally attached mesogens (17] [181, where the shape of the curve 
for 5 n had a semi-dome-like shape indicative of a power law form. More importantly, in the case consid¬ 
ered there, the decrease of 5 n and of the asphericity of the molecules was synchronous: both turned to 
zero at about T = 500 K. For the laterally attached mesogens considered here, we observe an essential 
delay in the decrease of two other order parameters, 5 h and Sq, as compared to 5 n. As one can see in 
figure]^ there is a narrow temperature range around T = 450 K, where the value of Sn dropped to about 
0.1 (typical of the isotropic phase), but 5 h and Sq are still almost the same as in the uColh morphology 
at r = 300 400 K. One may refer to this morphology as a weak hexagonal columnar (wColh) and the 
important feature of it is that while it lost the orientation order of mesogens, its hexagonal columnar 
structure is essentially preserved (see snapshots in the middle column of figure |^. One should note that 
the transition from uColh to wColh is gradual and the boundary between both shown in figurej^is rather 
for illustrative purpose. With further heating of the system, the values for the order parameters 5 h and 
Sc drop to their minima at approximately T ^ 480 K. The phase observed at the temperatures higher than 
that is a cubic phase (C), where the MPs adopt the form close to spherical due to an increased conforma¬ 
tional entropy. Its structure is shown in the right-hand column of figureand is the same as for the case 
of terminally attached mesogens I17II181 . 

The faster decay of 5 n with an in¬ 
crease of the temperature and, as a re¬ 
sult, the evidence of a new wColh mor¬ 
phology can be explained by the fol¬ 
lowing considerations. Let us go back 
to figure (a) and (b) and consider the 
difference in the allowed mechanisms 
for the mesogens reorientation in both 
cases. These are shown schematically 
in the same figure, frames (c) and (d) 
for the cases of terminal and lateral at¬ 
tachment, respectively. For the case of 
a terminal attachment, if we assume 
that the grafting point does not drift 
much, the only way for a mesogen to 
change its orientation is to bend the 
spacer [see, figure (c)]. There will 
be the energy penalty associated with 
this change due to the pseudo-angle 
potential term which is introduced to 
model the stiffness of the alkyl chain 
at a meso-scale level. This mechanism 
of bending the chain also exists for the 
case of laterally attached mesogens [figure[^(d), top chain with respect to the central sphere]. However, 
there also exists another mechanism of rotating a mesogen around the bond which connects it to the 
end monomer of a spacer [figure(d), bottom chains with respect to the central sphere]. This rotation 
costs no energy, as far as both bond length between a mesogen and the last monomer and the angle be¬ 
tween the orientation of a mesogen and the same bond stay unchanged. As a result, the rotational motion 
of mesogens in the mesogen-rich regions is more decoupled from the elastic energy of the spacers in 
the case of laterally attached mesogens compared to the case of their terminal attachment. This is, most 
likely, the main reason for a faster decrease of the nematic order with the raise of the temperature in 
the former case. We should note that the conformational changes in real systems involve more factors 
to be considered, as compared to the current coarse-grained modelling. In particular, both reorientation 
mechanisms depicted in figure[^(d) are likely to be realized via soft interactions such as torsional angles 
changes. Therefore, a more detailed analysis is needed to compare the energy penalty for both mecha¬ 
nisms of the mesogens reorientation discussed here. 

To conclude, we have shown the possibility of the formation of two novel morphologies, uColh and 


uCol h wCol h C 


Y 

A 

X 



Figure 5. (Color online) Shapshots showing the structure of 
different morphologies. Columns from left to right: uniaxial 
hexagonal columnar at T = 300 K (uColh), weak hexagonal 
columnar at T = 450 K (wColh) and cubic at T = 520 K (C) mor¬ 
phologies, all at the pressure of 100 atm. Rows from top to bot¬ 
tom: top view and side view for all respective cases (Cartesian 
axes are indicated in the figure). 



44001-5 












A.Y. Slyusarchuk, J.M. Ilnytskyi 


wColh for the MPs with lateral attachment of mesogens. The former is uniaxial in both the symmetry of 
MPs packing and in the arrangement of its mesogens. The latter is characterized by the same symmetry 
of MPs packing, while its LC subsystem is in the isotropic phase. Both morphologies differ in their LC 
properties from their counterpart for the MPs with terminal attachment of mesogens and are candidates 
for specific opto-mechanical applications. 
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HOBI M0p(|)0/10rM MeTdHdCTI/IHOK AeKOpOBdHI/IX I1J/15IXOM GiHHoro 
npi/ieAHdHH5i piAKOKpi/iCTd/iiHHi/ix rpyni KOMniorepHa ci/iMy/i5i 14151 

A.fO. CfiKDcapMyi^, 51.M. I/ibHi/mbKi/ibP 

^ HaqioHa/ibHUM yHiBepci/uer “/IbBiBCbKa no/iirexHiKa", By/i. C. BaH^epi/i, 12, 79013 /IbBiB, YKpaiHa 
^ iHCTi/uyT (^\3]AK]A KOH^eHCOBaHnx ci/icieM HAH yKpaiHH, By;i. I. CBCHqiqbKoro, 1,79011 /IbBiB, YKpaiHa 

PoBMfiHyTo MesocKoniMHy Mo^e/ib HaHopo3MipHi/ix MeraHacTi/iHOK (MM), 9\k\ cK/ia^aeibCfi is qeHTpa/ibHoi c(|)e- 
pi/i, AeKopoBaHoi no/iiMepHi/iMi/i /laHqwxKaMi/i is OImhi/im npi/ieAHaHH?iM c(t)epoi4i/i/iiHApiB. OcraHHi oni/icyK)Tb 
MesoreHHi (Hanp. i4i/iaHo6i(t)eHi/ibHi) rpyni/i. Bi/iKOHaHi CHMy/ifiqii 100 MM 3a ^onoMoroK) MO/iexy/uipHoi Ai^na- 
MiKH BKa3yK)Tb Ha icHyBaHH?! ABOX HOBi/ix Mop(t)o/ioriCi: uColh (i3 reKcaroHa/ibHi/iM croBnqeBi/iM BnaxyBaHH^iM 
MM ra ci/i/ibhi/im oahobIchi/im Bnop?iAKyBaHH?iM piAKOKpi/icra-niMHi/ix rpyn Ko/iineapno ao oci croBnqiB) ra wCol^ 
(i3 THM xe BHopjiAKyBaHHfi MM a^ne c/iaOKHM a6o BiAcyrniM piAKOKpi/icra/iiMHi/iM Bnop?iAKyBaHHJiM). Ko^nineap- 
HicTb Ai/^peKTopa ra oci CHMeipii croBnqeBoi (t)a3i/i y Mop(|)o/ioriT uCol^ BKa3ye na iT norenqinHO inmy ohti/iko- 
MexaniMHy peaKqiK) na 30BHiujHe 36ypeHH?i nopiBH^iHO i3 croBnqeBOK) (t)a30K) a^^ MOAe/ii i3 KinqeBi/iM npi/i- 
CAHaHHJiM piAKOKpHcra/iiMHi/ix rpyn. flonepeAHiCi ananis crpyKiyp o6ox (^as Bi/iKonano uj;i51xom AOC/iiAxenna 
napaMeipiB nop^AKy ra Bi3ya/ii3ai4ieK) po3TaujyBaHHri MM. BxaBano na po/ib piBHHX MexaHi3MiB peopieHiaqii 
Me3oreHiB 3a/iexHO BiA cnoco6y Tx npi/ieAHaHH?i. 

K/iiOHOBi c/ioea: HaHonacTi/iHKi/i, pifl,Ki Kpi/icra/ii/i, caMOBnopfi/^KyBaHHfi, Mo/ieKy/iBpna fli/inaMiKa 
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